Light-shift-induced quantum gates for ions in thermal motion 
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An effective interaction between trapped ions in thermal motion can be generated by illuminating them simul- 
taneously with a single laser resonant with the ionic carrier frequency. The ac Stark-shift induces simultaneous 
'virtual' two-phonon transitions via several motional modes. Within a certain laser intensity range these transi- 
tions can interfere constructively, resulting in a relatively fast, heating-resistant two-qubit logic gate. 
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Over the past few years, the implementation of a practi- 
cal quantum information processor has become a major goal 
for experimentalists across a wide range of disciplines |jl|]. 
Among the many candidate scenarios for attaining this goal, 
one of the most promising and well-known is the system of 
laser-cooled trapped ions, first suggested by Cirac and Zoller 
(CZ) . Several key features of their seminal proposal have 
already been experimentally demonstrated Meanwhile, 
the inherent difficulty of these experiments has stimulated the 
development of an array of alternative gate schemes [||-^ 
making increasingly ingenious use of the ion trap's physics. 

Here we combine the best features of two of these propos- 
als with an extra twist, constructing a two-qubit gate mecha- 
nism which is at once (i) robust against environmental heat- 
ing, (ii) relatively fast, and (iii) requires relatively few experi- 
mental resources, namely a single pulse by a single laser and 
no ground-state-cooling. Our first ingredient is the idea, due 
to S0rensen and M0lmer [^, of coupling different ions us- 
ing Raman-like two-phonon exchanges via a "data bus" mo- 
tional mode. Since the mode is only excited 'virtually', this 
allows the realisation of gates that are relatively insensitive to 
the ionic vibrational state, and can be implemented even in 
the presence of moderate motional heating This method's 
largest drawback [|[| is the very low switching rate of the 
gates, which are substantially slower than those obtained with 
CZ's method. This is undesireable, since it increases the quan- 
tum register's vulnerability to other sources of decoherence 
such as technical noise in the experimental apparatus. The 
problem has been tackled to an extent in ref. [01, where an 
elegant means of considerably increasing the gate speed is 
provided. However, this solution requires heavily entangling 
the internal and motional variables during the gate operation, 
making it more sensitive to heating. 

Our second ingredient is a method we have recently pro- 
posed for obtaining relatively fast gates between ground-state- 
cooled ions [^. The method reUes on the ac Stark-shift (light- 
shift) induced by driving an ion with a resonant laser. Nor- 
mally, this driving merely leads to Rabi flopping between the 
ionic states, realising a one-qubit gate. However, when the in- 
duced level-splitting is equivalent to exactly one motional en- 
ergy quantum, the influence of off -resonant motion-affecting 
transitions can gradually build up, leading to the exchange of 
excitations between the ion's internal and vibrational states. 
This can then be used to derive a fast CZ-hke two-qubit gate 



between any two ions in a linear chain. As with CZ's method, 
however, this scheme requires ground-state cooling and is 
highly sensitive to motional heating. 

In what follows we show a way of combining these ideas 
and their resulting benefits. Specifically, we propose to use the 
light-shift effect to drive 'virtual' two-phonon transitions via 
several motional modes simultaneously. Each mode functions 
as a parallel "data bus", with an overall gate resulting from 
the interference of various 'bus' paths. The gates have speeds 
intermediate between those in refs. [||] and [^, but remain in- 
sensitive to heating at all times. Moreover, only a single laser 
beam is needed, instead of the bichromatic illumination used 
in 1^,01 . As in the same beam can also be used to drive 
single-qubit gates. 
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FIG. 1 . A laser resonant with the carrier transition illuminates two 
ions, splitting each dressed level pair |±) by Ml. When this splitting 
is detuned from all vibrational frequencies Vq, states |H — ) 
and I — h) i {wq }) are degenerate with each other and non-degenerate 
with all other levels (for clarity, only levels of a single mode p are 
shown). Coherent oscillations are then induced between them, at a 
frequency independent of {uq}. This is due to the existence, for 
each mode, of 'virtual' two-phonon transition paths connecting the 
two states, via levels |++) |np — 1) and | ) \ np + 1). 

Let us begin with an intuitive description of our essential 
idea (Fig. Consider a chain of N ions, two of which are si- 
multaneously and equally illuminated by a laser resonant with 
the carrier frequency. At first sight, this simply causes these 
ions' internal states to flop periodically at the Rabi frequency 
il. Equivalently, each ion suffers a level-splitting of its (semi- 
classical) dressed states |±) — l/^/2{\g) ± |e)) by Kl. The 
presence of the trap, however, opens the possibility of tran- 
sitions between these dressed levels, via the absorption and 
emission of phonons Which transitions are energetically 
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allowed depends on the magnitude of the level-splitting. For 
example, when ^l is exactly equal to the frequency i/p of one 
of the ions' collective motional modes, then the manifolds 

{| ) K + 1) , K) , K) , \np - 1)} are 

degenerate and one expects transitions altering the phonon 
number in that mode. This is interesting, but inadequate for 
obtaining a motion-independent ion-ion gate. If, however, 
is slightly detuned from this value, then transitions involving 
changes in rip become off-resonant and are suppressed. One is 
left only with oscillations between | — h) \np) and |H — ) |np), 
which can be interpreted as occurring via the Raman-like ab- 
sorption and emission of 'virtual' phonons. 

Now, just as in ref. [^, each of these virtual transitions can 
occur via two different paths: one 'via' level |++) \np — 1), 
with amplitude proportional to rip and positive detuning Ap, 
and one 'via' level | ) |?ip + 1), with amplitude propor- 
tional to np + 1 and negative detuning — Ap. The overall 
transition amplitude will thus be the difference of these terms; 
remarkably, the ?ip-dependence should cancel out, and we can 
expect motion-state-independent transitions between the in- 
ternal states I — V) and |H — )! Of course, the ion system con- 
tains in fact several motional modes, each of which will lead 
to a separate resonance in the Rabi frequency. However, if 17 
is sufficiently detuned from all these resonances, the argument 
above holds simultaneously for all modes. In other words, the 
overall transition amplitude between | — h) and |H — ) should 
result from the interference of multiple Raman-like paths, two 
for each separate mode. Thus, the totality of modes are used 
together as a 'collective data bus' connecting the two ions. 

To see that all this is really the case, consider the Hamilto- 
nian of this system. Within the Lamb-Dicke limit 
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(where r\p is the Lamb-Dicke (LD) parameter of the p*'* col- 
lective mode of the ions and fip is the average number of 
phonons in that mode), this can be written as [0] 
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Hamiltonian becomes then 



and in our convention |e) 



(o) ■ The 
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+ e'^^'a],] + h.c. , 



(4) 



where we define Jpj. = Vjp^'j± ; Ap = — z/p, 

7p = + f p and primes indicate operators defined in the 
dressed-state picture. One can see that, when Ap — 0, the 
dynamics is dominated by terms of the form a'^^ap, which in- 
duce resonant collective excitations of the ions accompanied 
by the exchange of motional quanta. 

In what follows we are interested in the complementary 
regime far from these resonance conditions. Specifically, let 
us consider the limit where the detunings Ap of the Rabi 
frequency are large with respect to the resulting secular fre- 
quency of the time evolution, i.e. 



Ap > rijpfl = r]jp [vp + Ap) . 



(5) 



In this case, standard time-averaging arguments (see e.g. the 
Appendix in |^]) can be used to obtain a time-independent ef- 
fective Hamiltonian, given by 
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[ Jp_a|, Jp^Qp'j 



[jp_ap, Jp^a|] 
7p 



(6) 



Evaluating the commutators, we get after some algebra 

j;+ap] ^B'p-A'p + fci; [j^.Op, J^+at] ^ B'p + 



Ap + k2, where 



A'p - vipmp [\e'g') {g'e'l + \g'e') {e'g' 

r 

+ ivl-vl) [\g'e') {g'e'l -\e'g'){e'g'\ 



^v + 7i]Mp + v%) [\g'g') {g'g'\ 



(7a) 



(7b) 



+ h.c. (2) and where /ci, k2 are constants that may be disregarded by a 
suitable redefinition of the zero of energy. Thus 



Here f2 is the effective Rabi frequency of the laser-ion inter- 
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action (assumed to be equal for both ions), (pj is the laser's H' w —uo \\e' g') {g'e'\ + \g'e') (e'g'j] + 



phase at the position of the j*-^ ion, Vp is the frequency of 
the p*'* collective mode (p = 1 for the centre-of-mass (CM) 
mode, p = 2 for the 'breathing' mode, etc.) and rjjp is an 'ef- 
fective' LD parameter incorporating the relative displacement 
of the j*^ ion in the p*'' mode [p|,p^. We have also set h= 1. 
In what follows we will assume, with no loss of generality, 

01 - 02 = 

Following now a reasoning analogous to that of ref. [|^], we 
find that resonance conditions arise for particular values of il. 
This can be most easily seen by transforming into a dressed- 
state picture defined by the operator 



ptt ^i^' - -I) 



(8) 



where 
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V (t) = exp 



Ri (g) exp 



2z 



(3) 



Of course, this Hamiltonian holds only in the dressed pic- 
ture defined in Eq.(^. In the 'standard' picture (i.e., the 
one corresponding to Eq. (^), the time evolution operator 

is U{t) EE V\t) exp(-iti/^^^)y(0). 

In order to analyse this result, let us first consider the case 
= 2, where -qn = 7721 = \/ir]i2 = -v^?722 The 
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second term in Eq. ( [7b[ ) then vanishes, and H'^jj.- sphts into 
one term coupHng internal states \g'e') and \e'g'), and an- 
other affecting only \g'g') and |e'e'). Translating back into 
the "standard" picture, and including also the motional vari- 
ables, the first term leads to transitions between each pair 
of states |H — nin2) , | — h ?^l?^2), as expected from Fig. |l] 
(note that all time-dependent phases originating in V (t) can- 
cel out). Furthermore, as was also previously suggested, the 
frequency uj of the transitions is independent of the phonon 
numbers {np}. Finally, it is clear from Eqs. that this 

coupling results from the interference of contributions from 
both modes. The interference is constructive (both terms con- 
tribute positively to cu, leading to faster oscillations) when 
i^i < n < 1/2 = \f^v\ . Note also that in general the further 
is to Vq, the smaller is the amplitude in uj of the transition 
path via the corresponding mode. This effect might be useful 
to screen out 'bad' modes with relatively high decoherence 
rates, such as the centre-of-mass mode in multi-ion traps pl|l . 




Rabi frequency (n/ v) 
FIG. 2. Time required for the creation of a Bell state starting from 
the initial state |H — ), calculated from Eq. (^. Thin solid curves rep- 
resent the limits imposed by Eq. (^. The method should hold well 
in the region above both curves. 

Meanwhile, states |+ + n\n2) and | n\n2) acquire time- 
dependent phases, originating both from V (t) and from 
a'^j f- The latter do depend on the phonon numbers and can 
lead to unwanted correlations between internal and motional 
variables, and thus to decoherence of the internal state. These 
results are analogous to those obtained in ref. [|[| for the case 
where each ion is separately illuminated by a different laser, 
one shghtly detuned from the first red sideband frequency and 
the other from the first blue one. As in that case, here it turns 
out that the motion-state-dependent component of the evolu- 
tion can be cancelled by altering the Hamiltonian in the course 
of the gate operation itself (see below). 

As a result, maximally entangled states of the two ions can 
be produced after a time t\ — |7r/4w|. To see this, note that 
the evolution operator U{t\) is essentially a '\fSWAP^ gate 
with respect to the |±) basis: it turns the disentangled states 
|H — ) , I — h) into the maximally entangled (or 'Bell') states 

|/3±) EE (I+-) ± 1 1-+)) /\/2, while leaving |++) and | ) 

unaffected (up to motion-state-independent phases <^±(t]) 



due to V'^if)'). We stress again that this evolution occurs re- 
gardless of the initial motional state, and even regardless of 
whether it might be changing (due to heating effects) during 
the gate operation itself. Note further that, after a pulse of 
duration 2ti, |±=f) ^ i |t±), |±±) ^ exp(±if]Ti) |±±), 
and so all four standard basis states \g(e)g{e)) become max- 
imally entangled. In other words, a gate locally equivalent to 
a CNOT is realised. Our method therefore allows the imple- 
mentation of universal quantum logic on the ion chain, even 
in the presence of heating. The speed with which this is ac- 
complished can be seen in Fig. ^ where we use Eq. (^ to 
plot T\ as a function of 17, r\\\. The thin solid lines indicate 
the validity limits imposed by choosing a factor of 10 in Eq. 
(|^). We can conclude that the method should allow ri to be as 
low as a few hundred trap periods, a performance intermedi- 
ate between that in [pfl and that of the enhanced but heating- 
sensitive method in [0]. A significant difference with regard 
to these proposals is that our method should work in a dif- 
ferent parameter range, including a comparatively high Rabi 
frequency and a comparatively small LD parameter This may 
give it a substantial speed advantage in situations where ?;ii is 
forcibly small, such as a tight trap containing many ions. 

The trick which cancels the motion-state-dependent part 
of the evolution is to use a 'photon-echo' -like procedure [^], 
whereby midway along the time evolution one inverts the sign 
of the motion-dependent term in the Hamiltonian (in our case 
Bp in Eq.(|^)). As a result, phases acquired in the first half of 
the evolution are cancelled out exactly during the second one. 
In the present setup, this idea can be implemented in an ex- 
perimentally simple way, by suddenly shifting the laser phase 
by TT (using, for instance, an electro-optic modulator). To see 
this, let us return to eq.(^. Shifting (pj+n is equivalent 
to changing the sign of fi, which corresponds to exchanging 
the signs of the light-shifts suffered by |+) and |— )). A devel- 
opment analogous to eqs. (^-(||) results then in an effective 
Hamiltonian identical to i?eyy but with |e') and \g') every- 
where interchanged. It can be easily seen that this leaves the 
motion-independent term Ap intact, while changing the sign 
of Bp as required. In fact, following again ref. [||], in order 
to suppress motional heating it is desireable to perform the 
sign inversion several times during the course of the system's 
evolution (specifically, at frequency F — M/r, where r is the 
total time during which the laser is applied and M 1 is an 
integer). If F ^ oj, F, where T is the typical heating rate, 
then the effects of the motion-dependent terms are cancelled 
before heating can affect them [pj|]. 

This robustness is illustrated in Fig. |^, where we plot time 
evolution curves simulated using the full trap Hamiltonian 
(i.e., including all orders of the LD parameter) and allow- 
ing also for heating in the form of quantum jumps, described 
by jump operators ^Tpfipa and ^Tp{np + l)a^ [p4[]. (The 
heating rate Fi for the CM mode is assumed to be an order 
of magnitude greater than r2 [pi]]). The curves are averages 
obtained over 25 Monte Carlo runs, each containing an av- 
erage of 18.4 jumps. In the bottom curve, the ions' internal 
state is initially |i/'o) = |H — ), and we plot its squared over- 
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lap over time with the Bell state |/9_). A fidelity of 98% is 
achieved at the time j/iTi ~ 515 expected from Fig. ^. In 

the top curve, \tpo) = (|++) + | )) /\/2, and we plot the 

squared overlap with V'^{t)V{0) \ipo)^ the expected evolution 
in the absence of the motion-state-dependent term in Eq. (|^). 
The small deviation from 1 shows that this term has been ef- 
fectively suppressed by the switching of the laser phase. 
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FIG. 3. Time evolution of ions in thermal motion. In the bot- 
tom curve the initial internal state is l-i/io) ~ |H — ), and we plot the 
squared overlap over time with the Bell state In the top curve, 

\tpo) = (I++) + I )) /\/2, and we plot the squared overlap over 

time with yt(f)V'(0) I^Ao)- In both cases, the motional modes are 
initially in thermal states with mean phonon numbers ni = 1, n2 = 
0.1, undergoing heating at rates Fi — 10~^i^i, F2 = lO^^i^i. Other 
parameter values are 7711 — 0.025, Q, = 1.5i^i, F — i/i/SO. 

At first sight, it may appear that ever greater suppression of 
the motion-dependent evolution can be achieved by increas- 
ing the switching rate F. One must be careful, however, not 
to invalidate the approximations leading to Eq. (^. For in- 
stance, F must be limited by the Rabi frequency fl, in order to 
allow the rotating-wave-type averaging implicit in the deriva- 
tion of Eq. (||) to hold within each interval between phase in- 
versions. In fact, our numerical simulations indicate that this 
averaging only occurs fully for particular "resonant" values 
of F. For most other choices, the time evolution of the sys- 
tem gradually becomes degraded. Determining analytically 
these 'good' values is an open question; in practice, one would 
likely "tune" F until an appropriate value was found. 

In multi-ion chains with N > 2 our conclusions still hold, 
but high-frequency phase shifts become a necessity even in the 
absence of external heating. To see this, consider first that the 
LD parameters rjjp will generally differ from ion to ion [||], 
and so [Ap, Bp] ^ 0. In this case, inverting the sign of Bp at a 
low frequency (e.g., comparable to uj) is not sufficient to undo 
its effect on states |H — ) and | — h). However, if the frequency 
is much larger than uj then this cancellation does occur, since 
for very short timescales the effect of the commutator above 
becomes negligible [p^. As a result, once again Bp can be 
neglected, and the Hamiltonian reduces effectively to the first 
term in Eq. (^. A further generalisation to the case where 
m > 2 ions are illuminated together is also possible, as is the 
counter-intuitive possibility of coupling different pairs of ions 
simultaneously (cf the concluding remark in 1^). The latter 



case is achievable by illuminating each pair with a separate 
beam of sufficiently different intensity. In this case, following 
again the argument in Eqs. (!])-(§) shows that the Hamilto- 
nian effectively decouples into a sum of terms like in Eq. (|^), 
one for each pair. Thus, adding a further mind-boggling twist 
to the situation in [^], every single vibrational mode can be 
used collectively and simultaneously for different tasks, de- 
spite only ever being 'virtually' excited! 
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